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ABSTRACT 

In this research, thermal residual stresses were incorporated in an analysis of fiber-bridged 
matrix cracks in unidirectional and cross-ply titanium matrix composites (TMC) containing 
center holes or center notches. Two TMC were investigated, namely, SCS-6/Ti-15-3 and 
SCS-6/Timetal-21S laminates. Experimentally, matrix crack initiation and growth were 
monitored during tension-tension fatigue tests conducted at room temperature and at an 
elevated temperature of 200°C. Analytically, thermal residual stresses were included in a 
Fiber bridging (FB) model. The local R-ratio and stress-intensity factor in the matrix due to 
thermal and mechanical loadings were calculated and used to evaluate the matrix crack growth 
behavior in the two materials studied. The frictional shear stress term, r, assumed in this 
model was used as a curve-fitting parameter to matrix crack growth data. The scatter band in 
the values of r used to fit the matrix crack growth data was significantly reduced when 
thermal residual stresses were included in the fiber bridging analysis. For a given material 
system, lay-up and temperature, a single value of r was sufficient to analyze the crack growth 
data. It was revealed in this study that thermal residual stresses are an important factor 
overlooked in the original FB models. 
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Current crack length and unbridged initial crack length, m 
Superscript referring to crack centerline 
Cross sectional area of fiber and matrix, m 2 

Superscript referring to end line of slip region 

Hole diameter, m 

Crack growth coefficient constants 

Fiber and matrix modulus, MPa 


Composite longitudinal, transverse, and shear modulus, MPa 
Boundary correction factor for applied stress-intensity factor 
Correction factor for parallel cracks 
Weight function 

Applied stress-intensity factor range, MPa's/m 
Discrete stress-intensity factor range in matrix, MPa>/m 
Continuum stress-intensity factor range in composite, MPa*\/m 


Slip length, m 

Slip length of frictional shear stress reversal, m 
Distance between parallel cracks, m 
Superscript referring to maximum applied load 
Superscript referring to minimum applied load 
Range in bridging pressure, MPa 
Fiber radius, m 
Hole radius, m 


Composite stress ratio = 


Matrix stress ratio = 


/ ° 


max 

m 


Applied stress range, MPa 
Change in temperature, °C 
Fiber and matrix volume fractions 


Width of specimen, m 

Integration variable along crack from center, m 

Coefficient of thermal expansion for fiber and matrix, mm/mm/°C 


Continuum crack opening displacement range in composite, m 
Discrete crack opening displacement range in matrix, m 


Composite fracture surface energy 
Matrix fracture surface energy = 



MN/m 


MN/m 


Interfacial frictional shear stress, MPa 
Composite Poisson's ratio 
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0 f ,a m Axial stress in fiber and matrix, MPa 

o 1 l f* rm Axial thermal residual stresses in the matrix, MPa 

INTRODUCTION 

In several studies [1-4], fiber-bridged matrix cracking has been identified as a major fatigue 
failure mechanism in titanium matrix composites (TMC). Reductions of up to 38% in the 
composite longitudinal stiffness and static notched strength of cross-ply SCS-6/Ti-15-3 
laminates containing center holes have been reported due to fiber bridged matrix cracking [4]. 
Efforts to analyze the growth of these cracks using fiber bridging (FB) models have been 
undertaken in many investigations [5-10]. In these FB models, it is assumed that fiber-matrix 
debonding occurs as the matrix crack progresses past a fiber and a constant, but unknown, 
shear stress r acts along the debond length. The FB models combine a continuum fracture 
mechanics analysis and a micromechanics analysis to obtained stress-intensity factor solutions 
for fiber-bridged matrix cracks. 

In the research reported in [6], several of the original FB models were used to evaluate 
fiber-bridged matrix crack growth in unidirectional TMC. The frictional shear stress r 
assumed in the FB models was used as a fitting parameter for various data sets. In this study, 
it was shown that the FB models are an efficient and relatively simple engineering tool to 
conduct parametric analyses using the composite microstructural variables and may provide the 
framework for a damage tolerance methodology. There are, however, several shortcomings of 
the FB models that were pointed out in this study. 

Many aspects of fiber bridging crack growth were overlooked in the original FB 
models. The micromechanics analysis is based on a simplified one-dimensional force balance 
analysis of a single fiber in the wake of the matrix crack. Thus, these FB models do not 
account for the complex crack-tip singularity fields and the three-dimensional effects of crack 
front bowing around the fibers which can reduce the stress-intensity factor [11]. In addition, 
the assumption that a constant frictional shear stress acts along the fiber-matrix debond region 
is a simple engineering approximation. Variations in fiber-matrix debonding mechanisms were 
observed along the debond lengths in unidirectional TMC [6] which would yield different 
values in the shear stress. The interface stresses would also be affected by the wear of the 
debonded surfaces due to fatigue loading. Moreover, the large variation in the values of t 
used to evaluate the data in [6] indicated that a fundamental mechanism may be missing in the 
formulation of the FB models. The value of t required to fit various data sets depended on 
several factors, such as the applied stress and the crack extension length, and was not a 
material property [6], Consequently, the dependency of r on so many factors poses 
restrictions on using the FB models in a damage tolerance methodology. 

It was postulated that if the major damage growth mechanisms and the mechanics 
governing damage were considered in the FB model, r should represent a material constant. 
One important effect overlooked in the original formulations of the FB models is thermal 
residual stresses. Due to the mismatch in the mechanical and thermal properties of the 
constituents in TMC, thermal residual stresses are induced in the fiber and matrix during the 
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composite fabrication process and can effect the local and global fracture and fatigue behavior. 
Therefore, it is essential that thermal residual stresses be considered in the analysis of TMC. 

The objective of this research was to incorporate thermal residual stresses in the 
analysis of fiber-bridged matrix cracks in unidirectional and cross-ply TMC containing center 
holes or center notches. Two TMC were investigated, namely, SCS-6/Ti-15-3 and SCS- 
6/Timetal-21S laminates. Experimentally, constant amplitude, tension-tension fatigue tests 
were conducted under load control mode at a frequency of 10 Hz at room temperature and at 
an elevated temperature of 200 °C in order to determine the effects of thermal residual stresses 
on the crack growth behavior. Matrix crack initiation and growth were monitored using a high 
magnification closed circuit television (CCTV) system, a long focal length microscopic system 
with image acquisition capabilities, surface and edge replicas, optical microscopy, and 
scanning electron microscopy (SEM). Analytically, the matrix cracking observed in the two 
material systems was modeled using a fiber bridging (FB) model which was modified to 
include the effects of thermal residual stresses. The frictional shear stress term r assumed in 
this model was used as a curve fitting parameter to correlate matrix crack growth data. 
Experimental results were interpreted from the analysis made using the FB model. 


EXPERIMENTAL PROCEDURES 
Materials and Specimens 

Two TMC systems were studied in this research, namely, SCS-6/TM5-3 and SCS-6/Timetal- 
21S laminates. Both material systems were fabricated by Textron by hot isostatic pressing 
(HIPing) titanium foils between unidirectional tapes of silicon-carbide (SCS-6) fibers having a 
diameter of 0.14 mm. In the SCS-6/Ti-15-3 and SCS-6/Timetal-21S laminates, fibers were 
held in place with molybdenum and titanium-niobium cross weave wires, respectively. The 
composition of matrix materials used in the SCS-6/TM5-3 and SCS-6/Timetal-21S laminates, 
were Ti-15V-3Cr-3Al-3Sn and Ti-15Mo-3Al-2.7Nb-0.2Si, respectively. For the SCS-6/Ti-15- 
3 laminates, the lay-ups studied were [0]g and [0/90] s with a range of measured fiber volume 

fractions v f from 0.33 to 0.36. For the SCS-6/Timetal-21S laminates, the lay-ups studied 
were [0/90] s and [0/90] 2s with a range of measured fiber volume fractions v f from 0.36 to 

0.38. The material properties are listed in Table 1. Notice that some of the composite 
laminate properties were calculated using the VISCOPLY program developed by Bahei-El-Din 
[13]. 


Specimens were machined using a diamond wheel saw into straight-sided coupons with 
the 0° fibers in the loading direction. Each specimen was 152.5-mm in length. Center holes 
having diameters (d) of 6.3±0.3-mm were drilled ultrasonically in specimens having a width 
of 19.1-mm. Center slits having lengths (2a) of 5.08-mm, 7.62-mm and 8.89-mm were 
machined via electo-discharged machining (EDM) in specimens having a width of 25.4-mm. 
The SCS-6/Ti-15-3 specimens were tested in the as-received condition while the SCS- 
6/Timetal-21S specimens were heat treated at a temperature of 620°C for 8 hours in a vacuum 
and then fumaced cooled prior to testing. To enhance optical observations and surface 
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replication, the surface of each specimen was polished to obtain a flat and lustrous finish. 
Thin sheet aluminum end tabs having dimensions 19.1-mm x 19.1-mm were bonded on the 
ends of all specimens in order to prevent grip failure. Information on the specimens tested is 
listed in Table 2. 


Test Procedures 


Tension-tension fatigue tests were conducted under load control with R Conp 


= 0.1 



) 


at a frequency of 10 Hz using a closed-loop servo-hydraulic test machine equipped with 
hydraulic grips. Most of the test results discussed in this research were generated in previous 
studies [6,7]. In this study, one fatigue test was conducted at room or ambient temperature 
and another at an elevated temperature of 200° C using a quartz-lamp heater. Matrix crack 
initiation and progression were monitored and recorded in real time using a closed-circuit 
television system (CCTV) and a long focal length microscopic system capable of 
magnifications up to 325 x. Testing was periodically interrupted when increments in crack 
extension were observed to take surface replicas and to examine the specimens surface near the 
damaged region using a scanning electron microscope (SEM) and an optical microscope. 


Analytical Procedures 

Fiber bridged matrix crack growth was analyzed using the approach outlined in Figure 1. As 
shown, crack growth behavior is evaluated according to a version of the power law function 
introduced by Walker [12] which includes R-ratio effects: 


da 

dN 





AK ni 

p tO - m) 
K m/ J 


( 1 ) 


where C, n, and m are material constants for the neat matrix material, R m is the local stress 

min / 

ratio in the matrix ( R m = m / ) and A K m is the Mode I stress intensity factor range in the 

matrix. Thermal residual stresses can be included in the analysis through the quantities R ^ and 
A/T^.in Equation (1). 

First the stress ratio R m is calculated. Next, an iterative procedure is used to solve for 
AAT m with a FB model. In this procedure, a value of the frictional shear stress r is assumed. 
Equation (1) is then integrated numerically using Simpson's rule to obtain the crack length as a 
function of the number of fatigue cycles and compared with experiments. If the calculations 
and the experimental values do not agree, A K m is recalculated using another value of r . This 
iteration is continued until a value of t is found so that A K m yields an acceptable comparison 
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between experiments and calculations. Details of the calculation of R„ and A K m are provided 
next. 


Calculation of R m 

nun / 

The stress ratio in the matrix (R m = m / max ) was calculated using the VISCOPLY code [13], 

/ ® « 

a micromechanics based model. The VISCOPLY code combines the vanishing fiber diameter 
(VFD) model, a thermo-viscoplastic theory, and laminated plate theory to evaluate the rate- 
dependent, high temperature response of symmetric laminates. Using the VISCOPLY code, 
the average local stresses and strains in the constituents can be calculated and the global 
composite response can be predicted under complex thermomechanical conditions including the 
composite fabrication process. A more detailed description of the VISCOPLY code is 
provided in [13]. 

The thermo- viscoplastic properties of the constituents used in the calculation of R„ are 
reported in [13] for SCS-6/TM5-3 laminates and in [14] for SCS-6/Timetal-21S laminates. In 
order to include thermal residual stresses in the calculations, the composite fabrication process 
was analytically simulated by modeling the cooling down from an assumed stress-free 
temperature to room temperature A stress-free temperature of 555 °C was used for the SCS- 
6/Ti- 1 5-3 laminates and 620°C was used for the SCS-6/Timetal-21S laminates. For SCS-6/Ti- 
15-3 laminates, it was assumed that the stress-free condition occurs at one half the melting 
temperature of the titanium matrix, that is, at 555°C. For the SCS-6/Timetal-21S laminates, it 
was assumed that the residual stresses were removed during the heat treatment cycle at a 
temperature of 620°C. Values of the thermal residual stresses in the matrix in the 0° plies 
induced during cool-down of the composite fabrication a 7 *™ are listed in Table 1. 

After the cool-down, a global mechanical load cycle from S mm to was applied, 
and the corresponding stresses in the matrix in the 0° plies were calculated. Using these 
values of stress, the value of R m in the matrix in the 0° plies was calculated. Values of R m are 
listed in Table 2 for each specimen tested. As revealed by the results in Table 2, the values of 
7^, were considerably higher than the value of the applied R-ratio (Rc onp = 0.1). Thus, the 

value of Rc^p does not represent the actual stress ratio in the matrix and the constituent 7^, 
values listed in Table 2 should be used in the analysis. 

Calculation of AK m 

This section describes how the stress intensity factor range in the matrix AK m was calculated 
using a FB model which includes thermal residual stresses. As illustrated in Figure 2, it is 
assumed in the FB model that fiber-matrix debonding takes place as a matrix crack progresses 
past the fibers. The intact fibers in the wake of the matrix crack are idealized as a crack 
bridging pressure, A/>. An unknown but constant frictional shear stress r is assumed to act 
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along the debond length of the bridging fibers. A continuum fracture mechanics analysis and a 
micromechanical analysis are combined using discrete-continuum relations to obtain stress 
intensity factor solutions for fiber bridging matrix cracks of arbitrary size. 

Figure 3 is a flowchart showing the procedure used to calculate AAT M using a 
micromechanics analysis and a continuum fracture mechanics analysis. The crack opening 
displacement in terms of the bridging pressure obtained from both analyses are combined to 
yield a single non-linear integral equation in terms of the bridging pressure. Once the bridging 
pressure is known, A K m can be determined. 


Micromechanics Analysis 

As described in Appendix A, the crack opening displacement of the matrix, A5 m , can be 
expressed in terms the bridging pressure, A/? as: 

A5 m = XA p 2 (2) 


where: 


4 rv 2 f E f (E f v f +E m v m ) 

The terms in this equation are the composite microstructural parameters as defined previously 
in the nomenclature list. The frictional shear stress term, t, in Equation (3) is a difficult 
quantity to directly measure and was therefore used in this study as a fitting parameter to crack 
growth data. 

Thermal residual stresses do not effect the range relationship between the crack opening 
displacement and the bridging pressure given by Equation (2). Thermal residual stresses, 
however, do effect the relationship at the maximum and minimum loads: 


C =2X(p-“) 1 + 4X / ."“,.,£ / A7<o / -c,J (4) 

min 
m 

The second term in Equation (4) is the crack opening displacement due to a temperature 
change, A 7', where a f and a m are the coefficients of thermal expansion for the fiber and 

matrix, respectively. Equations (4) and (5) indicate that, when a f < a m , the absolute crack 
opening displacements will increase when thermal cool-down is accounted for. 
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Continuum Analysis 

The crack opening displacement in the composite, A5 Comp , is: 


A5 comp 


_ 2 _ 

E' 


EK Comp G(X,a,R,W)da 


( 6 ) 


where £' is an orthotropic material constant given by the laminate properties [15]: 



and A K Comp is the continuum stress-intensity factor range in the composite in terms of the 
unknown bridging pressure tsp\ 


bX Comf = AS^F apf - 2 SpG(x,a,R,W)dx (8) 

In this equation, A S is the applied far-field stress, F app is the applied stress boundary 
correction factor, a 0 is the initial unbridged crack length (hole radius or notch length). In 
addition, the term G(x,a,R,W ) in Equations (6) and (8) is the weight function for the 
bridging stress-intensity factor range. Compendiums of boundary correction factors and 
weight functions, F app and G(x,a,R,lV), are provided in [16-18] for a variety of cracked 

specimen configurations. For the two cases studied in this research, namely symmetric cracks 
emanating from a center notch or a center hole, expressions for F app and G(x,a,R,W ) are 

provided in Appendix B. It should be pointed out that the boundary correction factors and the 
weight functions in [16-18] were derived for isotropic materials and it was assumed that these 
terms can be applied to the orthotropic materials studied in this research 


Discrete-Continuum Relations 

The micromechanics and continuum analyses are combined using the relations: 



(9) 

aaw(aO 

(10) 
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The functions /(a 5 Com/1 ) and f{AK Comp ) are the discrete-continuum relations. The three most 

commonly used discrete-continuum relations were used in this study and are referred to as the 
Marshall, Cox and Evans (MCE) relation [8], the McCartney (MC) relation [9], and the 
McMeekings and Evans (ME) relation [10]. All three relations are reviewed in Appendix C. 
The crack opening displacements from the micromechanics analysis, Equation (2), and from 
the continuum analysis, Equation (6), are combined using discrete-continuum relations, 
Equation (9), to yield a single nonlinear integral equation in terms of the unknown bridging 
pressure Ap and the unknown frictional shear stress r . For example, using the ME relation 
[10], Ab m = Ab Comp , the bridging pressure can be expressed as: 



■I 


ASyfmF app - 2 ApG(x,a,R,WJdx 


) 


G( x,a,R,W )da 


(ID 


Assuming a value for t , the bridging pressure Ap is solved using a Newton-Raphson iterative 
procedure and a Simpson-type algorithm for numerical integration. Once the bridging pressure 
is known, the matrix stress-intensity factor range A K m can be obtained using Equations (8) and 
(10). Using the ME relation [10], AK m = A K Comp , the matrix stress-intensity factor becomes: 


AAT. = A 


-2 


ApG( x,a,R,W ) dx 


( 12 ) 


Appendix C outlines the derivation of A K m using the other two discrete-continuum relations. 


Application of FB Model to Cross Ply Laminates 

To apply the FB model to cross-ply laminates, it was assumed the fibers in the 90° plies had 
no effect on the growth of fiber bridged matrix cracks. In the micromechanics analysis, the 
fibers in the 90° plies were assumed to have no effect on the crack growth. The 90° plies 
were therefore modelled as matrix material alone. This implies that the fiber volume fraction 
for the cross plie laminates become a half of the actual value. Therefore the values of v f and v 
in Equation (3) were adjusted accordingly. It was also assumed that the fiber-bridged matrix 
cracks did not influence the far-field properties. Thus, in the continuum analysis, the actual 
fiber volume fraction was used in calculating the laminate properties for the orthotropic 
correction factor E' in Equation (7). It should be pointed out that fiber-matrix debonding has 
been observed in the 90° plies in cross-ply TMC resulting in a reduced composite modulus 
[4]. This phenomena, however, was not accounted for in this research. 
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Modeling Parallel Cracks 


Crack patterns such as that shown in Figure 4 where four cracks initiate and progress from a 
center hole have been commonly observed in TMC specimens [4,7]. A finite element 
analysis conducted by John, et. al [19] found that the stress-intensity factor for fiber bridged 
matrix cracks in TMC (SCS-6/Ti-24- 1 1) were reduced by 20% when dual parallel cracks 
similar to the schematics in Figure 4 were modeled compared to a single crack. In order 
model this effect in the current model, AK m was multiplied by a correction factor for two 
parallel cracks [18]: 


/, u , = 1-0.293 4 {|-(1- 1 )*) 


(13) 


where: 


a 

s = 

a + h 


(14) 


In this equation, a is the length of the cracks and h is the distance between the cracks. The 
measured values of h for each specimen are listed in Table 2. The values of F par ranged from 

1 to which bound the value of 0.8 obtained by John, et. al [19] in their finite element 
analysis. 


RESULTS AND DISCUSSIONS 

The matrix crack growth behavior observed in the two TMC was analyzed using a FB model 
with thermal residual stresses. The frictional shear stress term t in the FB model was used as 
a curve-fitting parameter to matrix crack growth data reduced using the crack growth law 
defined by Equation (1). The value of t was determined for each of the three discrete- 
continuum relations: the MCE [8], MC [9] and ME [10] relations. Results are presented in 
the subsequent sections according to the lay-up and notch geometry. The effects of thermal 
residual stresses on the values of r used in fiber bridging analysis are discussed. 


Baseline Matrix Crack Growth Curves 

The matrix crack growth curves used in the analysis for the Ti-15-3 and Timetal-21S materials 
are shown in Figures 5 and 6, respectively, with the corresponding material constants defined 
in Equation (1) . The value of m in Equation (1) represents the amount of shift in the crack 
growth curve due to the R-ratio. In order to evaluate m, crack growth data must be available 
for more than one R-ratio. For the two matrices considered in this study, the crack growth 
data was available for only one R-ratio, R m = 0. 1 Therefore, the value of m was calculated 
using Ti-6-4 crack growth data available at several R-ratios (unpublished data, R. Everett, 
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NASA Langley Research Center, 1993). It was assumed that the effect of R-ratio on the crack 
growth behavior would be the same for the three titanium alloys; thus, the same value of m 
was used. 


Unidirectional Center Notch Specimens 

The crack length as a function of the number of cycles at several applied stress levels was 
calculated using the value of r which best fit the experimental results obtained in [6] for two 
center notched specimens of [0]g SCS-6/Ti-15-3 tested at room temperature, Figure 7. In both 

specimens, a single crack initiated and grew from each tip of the machined slit with no fiber 
breakage. For one specimen, the applied stress level was increased after a prescribed number 
of fatigue cycles, while for the other specimen, the applied stress was held constant. In the 
former case, the crack growth rate increased as the applied stress level increased. For both 
specimens, the crack growth rate decreased as the crack length increased due primarily to the 
observed fiber bridging phenomena. 

As shown in Figure 7, when thermal residual stresses were included in the analysis, 
good agreement was obtained between the experiments (open symbols) and calculations (filled 
symbols) made using the MCE relations [8] and a single value of the frictional shear stress, r 
= 20 MPa. These results are quite encouraging since the value of r required to correlate the 
crack growth data when thermal residual stresses were neglected was a function of the applied 
stress level [6]. Similar results were obtained when the crack growth data was reduced using 
the MC [9] and the ME discrete-continuum relations [10]. A value of r = 15 MPa was used 
for the MC relation and a value of r =60 MPa was used for the ME relation. For all three 
discrete-continuum relations considered, a single value of r was sufficient to correlate the data 
when thermal residual stresses were included. 

Using the same value of r used to reduce the data shown in Figure 7, the debond or 
slip length as a function of distance from the first intact fiber was calculated and compared in 
Figure 8 to experimental results obtained in [6] for a center notched [0] 8 SCS-6/Ti-15-3 

specimen. After fatigue loading, the surface of a specimen was polished to the mid plane of 
the outer ply of fibers. The mechanisms of fiber-matrix interfacial debonding were identified 
and the debond lengths were measured using the scanning electron microscope (SEM) [6]. 
Fiber-matrix debonding was assumed to exist where cracking between the carbon rich coatings 
and fracture of the brittle reaction zone was observed. The measured debond lengths were 
assumed to be equal to the entire damage process zone in the interface region. 

In the FB model, the slip length is defined as the region over which the frictional shear 
stress acts. In the present study, the slip length is assumed to be equal to the debond length, 
and the measured debond lengths are compared to calculated slip lengths. As shown in Figure 
8, good agreement was obtained between the measured and calculated debond length data 
using the value of t used to consolidate the crack growth data in Figure 7. The calculations 
shown in Figure 8 were made using the MCE relation. When thermal residual stresses were 
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neglected, there was no correlation in the values of r required to fit the debond length and 
crack growth data [6]. 


Cross-Ply Center Notch Specimen 

The crack length as a function of the number of cycles at room temperature and at an elevated 
temperature of 200°C were calculated and plotted in Figure 9 using a value for r which best 
fit the experimental results for two center notched specimens of [0/90] s SCS-6/Timetal-21S. 
In both cases, only one crack initiated and grew from each tip of the machined slit with no 
fiber breakage. As a result of the fiber bridging phenomena, the crack growth rate decreased 
as the crack length increased. The crack growth rate increased as the applied temperature 
increased. 

As shown in Figure 9, good agreement was obtained between experiments (open 
symbols) and the calculations (filled symbols) using the MCE relations and values of r =50 
MPa and 15 MPa for the room temperature and elevated temperature test results, respectively. 
For the elevated temperature test, the room temperature properties of the titanium matrix were 
used in the fiber bridging analysis. It was assumed that the mechanical properties of the 
matrix were the same at both temperatures since the modulus of Timetal-21S at 200°C was 
only 3% less than that at room temperature. Hence, the difference in the calculated crack 
growth behavior at the two temperatures is due primarily to the differences in thermal residual 
stress states which effects the calculated values of R m and AK m . Using the MC relation, r = 
50 MPa and r = 15 MPa for the room temperature test and elevated temperature test, 
respectively. Using the ME relation, t = 100 MPa and r = 32 MPa for the room 
temperature test and elevated temperature test, respectively. 

If the mechanical properties of the matrix are assumed to be identical at both 
temperatures, only R M and A K m affect the crack growth rate defined by Equation (1). At 
elevated temperatures, the axial thermal residual stress in the matrix o T ^ rm is lower than at 
room temperature which results in lower values of R m . For the room and elevated temperature 
tests, the calculated local R-ratios in the matrix were R m = 0.75 and R m = 0.64, respectively. 
From Equation (1), for a given value of A K mt the crack growth rate increases with an increase 
in the value of R m . Based solely on the values of R m for the two test conditions, the crack 
growth rate should have been higher for the room temperature test than the elevated 
temperature test. Results shown in Figure 9, however, show the opposite trend; the crack 
growth rate data was higher for the elevated temperature test compared to the room 
temperature test. 

Another factor affecting the crack growth rate is the stress-intensity factor range AK M . 
Thermal residual stresses affect the value of A K m . The lower value of r for the elevated 
temperature test translated into a higher value in the stress-intensity factor range as shown in 
Figure 10. At the elevated temperature, the value of A K m was higher than at room 
temperature. Using the crack growth law defined by Equation (1), for a given value of R m , 
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the crack growth rates increases with an increase in the value of AK m which agrees with the 
results shown in Figure 9. As a result of the lower thermal residual stresses at the elevated 
temperature test, the value of A K m increased while the value of R m decreased. The lower 
values of r for the elevated temperature results suggests that the radial thermal residual 
stresses acting across the debonded fiber-matrix interfaces in the 0° plies were lower. These 
radial stresses "choke" the fibers which enhances the shear transfer between the fiber and 
matrix. By lowering the radial thermal stresses, the amount of shear transfer and displacement 
constraint between the fiber and matrix were lowered thereby increasing the stress-intensity 
factor. The increase in the value of A K m had a more profound effect on the crack growth rate 
than the decrease in the value of R m . 

m 


Cross-Ply Center Hole Specimen 

The crack length as a function of the number of cycles at several applied stress levels was 
calculated using a value for r which best fit the experimental results obtained in [7] for room 
temperature tests of cross-ply SCS-6/Ti-15-3 and SCS-6/Timetal-21S laminates as shown in 
Figures 11 and 12, respectively. The calculations shown in Figures 11 and 12 were reduced 
using the MCE relation. As shown in these figures, when the thermal residual stresses were 
included in the analysis, good agreement was obtained between the calculations (filled 
symbols) and the experiments (open symbols) using a single value of the frictional shear stress 
for each laminate, r = 10 MPa and r = 50 MPa, respectively. Using the MC relation, the 
values of frictional shear stress needed to fit the data were r = 10 MPa and r = 50 MPa for 
the SCS-6/TM5-3 and SCS-6/Timetal-21S laminates, respectively. Using the ME relation, 
the values of frictional shear stress needed to fit the data were t = 25 MPa and r = 100 MPa 
for the SCS-6/Ti-15-3 and SCS-6/Timetal-21S, respectively. 


Effect of Thermal Residual Stress on the Frictional Shear Stress 

The effect of thermal residual stresses on fiber bridging analysis is shown by the values of the 
frictional shear stress term required to reduce the crack growth data, listed in Table 2. The 
values of r listed in this table were determined using the ME relations with and without 
thermal residual stress. For the unidirectional SCS-6/Ti-15-3 specimens tested in [6], the 
value of r required to fit crack growth data varied from 12 MPa to 41 MPa when thermal 
residual stresses are neglected (i.e. original fiber bridging analysis). However, when thermal 
residual stresses were included in the analysis, a single value of frictional shear stress, r = 60 
MPa, was able to fit the data. Thus, for a given material system, lay-up, and temperature, a 
single value of r was sufficient to analyze the crack growth data. Thermal residual stresses 
are an important factor not included in the original formulations of the FB models. 
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CONCLUDING REMARKS 


Thermal residual stresses were implemented into the analysis of fiber-bridged matrix crack 
growth in notched titanium matrix composites (TMC). The matrix crack growth behavior in 
unidirectional and cross-ply SCS-6/Ti-15-3 and SCS-6/Timetal-21S laminates containing center 
holes or center notches was analyzed. Results from the analysis revealed that tensile axial 
stresses in the matrix were induced during composite fabrication and heat treatment resulting in 
higher constituent R-ratios in the matrix (i.e., ratio of stresses at the minimum and maximum 
applied loads) compared to the continuum R-ratio based on the applied load (i.e., ratio of the 
minimum and maximum applied load). Thus, the continuum R-ratio of the applied load is not 
appropriate in the analysis of fiber-bridged matrix crack growth; the R-ratio based on the 
matrix stresses should be used since the crack growth does indeed take place in the matrix 
material.. Thermal residual stresses were incorporated in the formulation of a fiber bridging 
(FB) model which combines a continuum fracture mechanics approach and a micromechanics 
approach to yield stress intensity factor solutions for fiber bridged matrix cracks. The 
frictional shear stress term r assumed in the FB model was used as a fitting parameter to 
matrix crack growth data. 

Tension-tension fatigue tests were conducted at room temperature and at an elevated 
temperature of 200°C. For a given material system, and lay-up, and temperature, a single 
value of r was sufficient to reduce the crack growth data. The crack growth rates were higher 
for the elevated temperature test compared to the room temperature test. The fiber-bridging 
analysis indicated that the increase in crack growth rates was due to the lower thermal residual 
stresses which resulted in an increase in the stress-intensity factor range in the matrix. 
Incorporation of the residual thermal stresses showed their effects in two significant ways: (1) 
The stress ratio in the matrix increases by including the effects of residual stresses over that of 
the applied stress ratio. The residual stress in the matrix is a function of the fiber volume 
fraction and lay-up. Therefore these variations will be reflected in the matrix stress ratio. (2) 
The residual stresses change as a function of temperature. The axial stress ratio in the 0° ply 
matrix material goes down as would the radial residual stresses as the temperature is increased. 
The effect of the lowered radial stresses is to directly lower x because less frictional stress will 
be transferred between the matrix and fiber. When testing at elevated temperatures, it was 
found in this study that the lowering of the frictional shear stress x due to the decrease in radial 
residual stresses is more profound than the lowering of the stress ratio in the matrix material. 

Results from this study showed that thermal residual stresses are an important 
fundamental factor that was lacking in the original FB models. By including thermal residual 
stresses in the FB model, fiber bridging analysis more accurately models the fiber bridging 
phenomena observed in TMC. 
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Table 1. Material Properties: 
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Table 2. Specimen loading history and fiber bridging analysis parameters: 
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Frictional shear stress required to analyze data using ME relations [10] when thermal residual stresses are neglected 
Not applicable 

Hole diameter to width ratio (d/W) 

Frictional shear stress required to analyze data using ME relations [10] neglecting the effect of parallel cracks 
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Figure 1. Analysis of Fiber-bridged matrix cracks 
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Figure 2. Continuum fracture mechanics approximation of matrix crack growth process. Fibers bridging crack are idealized as a 
closure pressure. Governing equation is formulated by superposition of a crack subjected to a far-field applied stress 
and a crack subjected to a closure pressure. 




Figure 3. Calculation of Stress-Intensity Factor using Fiber Bridging Model 
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Figure 4. Schematics of 




Figure 5. Room temperature crack growth rate data for Ti-15-3. R-ratio = 0.1 




Figure 6. Room temperature crack growth rate data for TimetaI-21S. R-ratio = 0. 1 




SCS-6/TM5-3 Center Notch 
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Figure 7. Crack growth data for center notched [0]g SCS-6/Ti-15-3 laminates at several applied stress levels. 
Comparison between experiments and calculations made using Fiber Bridging Model. 



SCS-6/TM 5-3 Center Notch 
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Figure 8. Debond length data for a center notched [0] g SCS-6/Ti-15-3 specimen subjected to 50,000 cycles. 
Comparison between experiments and calculations made using Fiber Bridging Model. 



SCS-6/Timetal-21 S Center Notch 
[0/90] 8 2a/W = 0.2 

V* = 0.36 R. = 0.1 
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Figure 9. Crack growth data for center notched [0/90] s SCS-6/Timetal-21S laminates at two temperature levels. 
Comparison between experiments and calculations made using Fiber Bridging Model. 



SCS-6/Timetal-21S Center Notch 

[0/90], 2a/W = 0.2 
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laminates at two temperature levels. 



CO CO 
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Figure 11. Crack growth data for center hole [0/90] s SCS-6/Ti-15-3 laminates at several applied stress levels. 
Comparison between experiments and calculations made using Fiber Bridging Model. 






o 
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Figure 12. Crack growth data for center hole cross-ply SCS-6/Timetal-21S laminates at several applied stress levels. 
Comparison between experiments and calculations made using Fiber Bridging Model. 


APPENDIX A - MICROMECHANICS OF FIBER BRIDGING 


The relation between the crack opening displacement and the bridging pressure is 
derived in detail for each region of fatigue loading shown in Figure Ala. A free body diagram 
of a fiber in the wake of the matrix crack for each loading phase is shown in Figures Alb and 
Ale. A constant temperature change, A7', was applied prior to the mechanical loading. 


Loading: O-A 

During the initial loading along O-A (Figure Ala), the frictional shear stress, r, is 
constant in the slip region, f. , and opposes the fiber stress, Figure Alb. This free body 
diagram is constructed assuming the maximum load is just reached (Point A in Figure Ala). 
Along line BB' (end of slip region), the strains in the fiber and the matrix are equal: 




(Al) 


\m*x / ynnx 

a f ) and [a m ) are the fiber and matrix stresses along line BB' at the maximum 

applied load (Point A in Figure Ala). Overall equilibrium in the matrix in the slip region 
requires: 


urn 4.-2™*= 0 (A2) 

where A m is the cross-sectional area of the matrix. Equilibrium in the fiber in the slip region 
requires: 



A f + 2 TT/f 


(«r"P 




(A3) 


where (a^’P is the fiber stress along the matrix crack centerline AA’ and A f is the cross- 

seclional area of the fiber. Combining Equations (Al), (A2), and (A3), the fiber stress along 
the crack centerline A A' at the maximum load is: 



(A4) 


where: 
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(A5) 




1 + - 


E f u 


f y f 


The total extension of the matrix in the slip region due to the mechanical and thermal loading: 





2 Terry 

XX 


dy + a m A7’f 


rrrt 1 
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+ a m AH 


and of the fiber: 


(A6) 
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+ u 
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2rrry 

XX 


\dy + a f ATE - 



Trr (. 2 
+ 


E f A f 


+ a f ATl 


(A7) 


Combining Equations (Al), (A2), (A6), and (A7), the crack opening displacement at the 
maximum load, 5^“: 



2 


+M'f.(ot / -a„) 

l.jt 


(A8) 


Substituting Equation (A4) into (A8) yields the relation between the fiber stress along the crack 
centerline and the discrete crack opening displacement at the maximum applied load: 


h 


mux 

m 


1 



4tEj1) 



(A9) 


Unloading: A-B 

Upon unloading, the frictional shear stress reverses direction within a length, t R> in the 
slip region. The free body diagram, Figure Ale, is constructed assuming the minimum load 
is just reached (Point B in Figure Ala). The procedure to determine the relation between the 
fiber stress and discrete crack opening displacement is similar to that at maximum load. Along 
line BB' (end of slip region), the strains in the fiber and the matrix are equal: 
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E. 


E. 


(A10) 


Overall equilibrium in the matrix in the slip region requires: 


(o f )""" A m -2Trr{e-t H ) + 2TTrt H =0 


(All) 


and in the fiber: 


(u ) mn A f f 2itTr{i-^ H )-2rrr(. H = (<r A /') A f 


(A 12) 


Combining Equations (A 10), (All), and (A 12), the fiber stress along the crack centerline at 
the minimum load is: 




(A 13) 


Substituting Equation (A4) into (A 13): 




(A 14) 


where Ao f - (a , ) ~\o f j • The total extension of the matrix in the slip region at the 

minimum load is: 


J * f e '( ( 
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(A 15) 


which reduces to: 
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( „BB : \ nun 0 
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(A16) 


The extension of the fiber at the minimum load is: 
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(A 17) 


which reduces to: 


* nun 

U m ^ nun 
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/ ««• ) n,in « 


7'7 


(A 18) 


Combining Equations (A 10), (All), (A 16), and (A18): 




(A19) 


By substituting Equation (A8) into (A 19): 


rh f 


(A20) 


where A5 m = 5 '7 - 5™". Combining Equations (A14) and (A20) yields the relation between 
the change in fiber stress along the crack centerline and the change in discrete crack opening 
displacement from points A to B in Figure Ala: 


A5. 


-(a of ) 1 

4 tE, 7j 


(A21) 
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The unknown bridging pressure, Ap, and the fiber stress along the matrix crack 
centerline, A o A /\ are related through the fiber volume fraction as follows: 


A/> ~ if fAo A f A (A22) 

This equation is valid for fibers are closely spaced relative to the crack length. Using 
Equations (A21) and (A22), the matrix crack opening displacement is related to the unknown 
bridging pressure: 

AS m - XA p 2 (A23) 

where: 

X (- - (A24) 

Arvjhj T/ 


Reloading: B-A 

When reloading back to A, the constant frictional shear stress in reverses direction. 
Upon reaching A, the shear stress within the entire slip length, i , is in the same direction. 
The free body diagram in this case is the same as that shown for the original loading case, 
Figure Alb. Consequently, the discrete crack opening displacement when reloaded to Point A 
is identical to Equation (A9). Moreover, the change in displacement from Point B to A in 
Figure Ala is the same as that from Point A to B (Equation (A23)). 
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APPENDIX B - BOUNDARY CORRECTION FACTOR AND WEIGHT 

FUNCTION 


For two symmetric cracks progressing from a center hole of radius R in an isotropic 
material, the applied stress boundary correction factor is [16]: 


‘ 'n 


where: 


and: 


1 + 0.358— - + 1.4251 
a 


O' 



(Bl) 


(B2) 


K 



(B3) 


For two symmetric cracks emanating from a center notch, F app is obtained by substituting R = 
0 into Equations (B2) and (B3). 

The weight function, G(x,a,R,W), for two symmetric concentrated forces, P, applied 
to the cracks surface in a center hole specimen configuration [16]: 


G( x,a,R,W ) = GF 


(B4) 


In this equation, G cn is the weight function for a center notched specimen configuration 
subjected to a pair of forces, P, on the crack surface and F h is the correction factor to account 
for a center hole boundary. The weight function for a center notched specimen is given by 
[181: 



(B5) 


where: 
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The correction factor to account for the center hole of radius R is given by [17]: 
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APPENDIX C - DISCRETE-CONTINUUM RELATIONS AND 
GOVERNING EQUATIONS 


The continuum fracture mechanics analysis and the micromechanics analysis are 
combined using the following two discrete-continuum relations: 


M. = /(a 5 Cw ) 

(Cl) 

a*„=/(a O 

(C2) 


Differences among the original fiber bridging models [8-10] arise in defining these two 
discrete-continuum relations as discussed next. 


Crack Opening Displacements 

Marshall et al. [8] and McMeeking and Evans [10] assumed that the discrete and the 
continuum crack opening displacements were equal: 

A5 m =A h Comf (C3) 


McCartney [9] reported that A5 Camp is actually the change in displacement within the slip 

region (change in displacement between lines AA' and BB', Figures Alb and Ale) and should 
be related to A 8„ as follows: 


M Comp V 


where: 


V = 


1 + 


V/ ] 

E m V m ) 


(C4) 


(C5) 


Stress-Intensity Factor 

The relation between the matrix and composite stress-intensity factors ranges was 
established by Marshall et al. [8] using the ratio of the matrix modulus to composite modulus: 
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(C6) 


A*. = ^ A* r _ 

This equation assumes that the near-tip strains in the composite and in the matrix are 
compatible. In a different approach, McCartney [9] used an energy balance to relate AK m to 
AK Comp by assuming: 


7 Comp ^ ni7m 


(C7) 


K 2 / K 2 / 

where y Comp = Com yj? is the fracture surface energy of the composite and y m = is the 

fracture surface energy of the matrix. Rewriting Equation (C7) in terms of stress intensity 
factor ranges: 


AA_ --= 




(C8) 


Finally, it was suggested by McMeeking and Evans [10]: 


- AX, 


Comp 


(C9) 


In this paper, the discrete-continuum relations derived by Marshall, et al. [8] (Equations (C3) 
and (C6)), McCartney [9] (Equations (C4) and (C8)) and McMeeking and Evans [10] 
(Equations (C3) and (C9)) are referred to as the MCE, MC, and ME relations, respectively. 


Governing Equations 

The crack opening displacements from the micromechanics analysis, Equation (2), and 
from the continuum analysis, Equation (6), are combined using discrete-continuum relations, 
Equation (Cl) to yield a single nonlinear integral equation in terms of the unknown bridging 
pressure, Ap, and the unknown frictional shear stress r. Using Equation (C3), the MCE [8] 
or ME [10] relation: 
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I ASj^F a „- 2 ApG( x,a,R,W )dx G(x,a,R,W)da (CIO) 


Using the MC relation [9], namely, Equation (C4): 


AS^F app - 2 ApG(x,a,R,W)dx \G(x,a,R,W )da (Cll) 


Assuming a value for t, the bridging pressure, Ap, is solved using a Newton-Raphson 
iterative procedure and a Simpson type algorithm for numerical integration. Once the bridging 
pressure is known, the matrix stress intensity factor range, A K m , can be obtained using 
Equations (8) and (C2). Using the MCE relation [10], namely, Equation (C6): 


A K m - ASyfxaF - 2 ApG(x,a, R, W)Jx 

K 


(Cl 2) 


Using the MC relation [9], namely, Equation (C8): 


A K m = I— ~ AS-JiaF - 2 ApG(x,a,R,W)dx 
\ v J'i 


(Cl 3) 


Using the ME relation [10], namely, Equation (C9): 


AK m = ASjniF ' -2 I Aj*G(x,aJi,W)dx 


(C14) 
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Fiber Stress 
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Figure Al. Micromechanics of fiber bridging: (a) Fiber stress as a function of crack opening displacement; 
(b) Free body diagram of fiber in wake of crack at maximum applied load; and (c) Free body 
diagram of fiber in wake of crack at minimum applied load. 
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